On the electrostatic potential profile in biased molecular wires 
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The potential profile across a biased molecular junction is calculated within the framework of a 
simple Thomas-Fermi type screening model. In particular, the relationship between this profile and 
the lateral molecular cross section is examined. We find that a transition from a linear potential 
profile to a potential that drops mainly near the molecule-metal contacts occurs with increasing 
cross section width, in agreement with numerical quantum calculations. 
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I. INTRODUCTION 

Molecular wires are molecules that bridge between 
metallic leads to form a nano-conductor whose current- 
voltage characteristic reflects the electronic structure of 
the molecule. The conductance may be controlled via 
its dependence on molecular properties. Equally impor- 
tant is the use of such molecular junctions as capacitive 
elements in nano-electronics. 

Understanding the behavior of such junctions under 
potential bias is a prerequisite for elucidating their trans- 
port properties. The importance of this issue for the 
conductance of molecular junctions was recently empha- 
sized by Datta and coworkers jj], |2), who have shown, 
within a simple Extended-Huckel (EH) model for a,a'- 
xylyl dithiol bridging between two gold leads, that the 
potential profile (imposed on the molecule as input to 
the EH calculation) had a profound effect on quantitative 
as well as qualitative aspects of the calculated current- 
voltage characteristic. The best fit to experimental re- 
sults was obtained from a model that assumed (a) a flat 
potential profile in the interior of the molecular bridge, 
i.e. the potential drop occurs only at the molecule-lead 
contacts and (b) a symmetric distribution of the poten- 
tial drop at the two contacts, i.e. for a total voltage <& 
the potential drops at each molecular edge by $/2. 

This picture is supported by a recent model calcula- 
tion by Mujica et al. ||, where the Schrodinger equation 
(on the Hartree level) was solved in conjunction with the 
Poisson equation to yield both the electronic structure 
and the spatial distribution of the electrostatic potential 
H . It was found that beyond a screening distance of the 
order of 1-3 atomic lengths the potential is flat along the 
model molecular chain. 

Ab initio calculations with open system boundary con- 
ditions reveal a different picture: Lang and Avouris |^| 
have found for a chain of seven carbon atoms connect- 
ing between jellium leads that a substantial part of the 
voltage drop occurs along the carbon chain itself. Damle 
et al. H have obtained similar results for a chain of gold 
atoms as well as for a molecular wire — phenyl-dithiol 
bridging between two gold electrodes 0. In an earlier 
work, Pernas et al. Q have determined that the poten- 
tial along a model molecular wire is flat in the absence 
of scattering centers, however these authors have derived 



the local potential from a questionable local charge neu- 
trality condition. 

Recently, Weber et al. || have considered the volt- 
age profile across 9,10-Bis((2'-para-mercaptophenyl)- 
ethinyl)-anthracene coupled to two AU29 clusters. Their 
density functional theory calculations thus go beyond the 
assumption of a structureless metallic electrode and take 
into account the specific properties of the bond between 
the molecule and the gold atom in its local environment. 

On the experimental side, Bachtold et al. jl0| have 
used scanned probe microscopy of both single- walled and 
multi- walled carbon nanotubes (SWNT and MWNT, re- 
spectively) to measure the potential drop along such nan- 
otubes connecting between two gold electrodes. They 
find an approximately linear drop of the potential in a 
MWNT of diameter 9nm while for a metallic SWNT bun- 
dle of diameter 2.5nm the potential is flat beyond the 
screening regions at the tube edges. It should be empha- 
sized that these experiments cannot be related directly to 
the calculations discussed above. The nanotube length is 
a few microns and impurity and defect scattering may be 
effective as is most certainly the case in the MWNT mea- 
surement. The flat potential seen in the metallic SWNT 
measurement is in fact a remarkable observation imply- 
ing a very long mean free path (> 1/xm) for electrons in 
these room temperature structures. 

It is clear from the studies described above that while 
the computational methodology for evaluating the poten- 
tial distribution on a biased molecular wire is available, 
a complete understanding of the way this distribution 
behaves in different wires is lacking. In this respect sim- 
ple models that focus on generic properties of conducting 
constrictions including molecular wires are useful. The 
calculations of Pernas et al. ^ provide such a model that 
is however hampered, as already stated by the restriction 
of local charge neutrality. The calculation of Mujica et al. 
H is also based on a generic molecular model, however, 
by using a 1-dimensional Poisson equation for the electro- 
static potential these authors tacitly assume a molecular 
bridge whose lateral dimension is much larger than the 
screening length. In view of the fact that the width of 
molecular wires is often just a few angstroms, such an 
assumption is overly restrictive. Clearly, the magnitudes 
of the lateral width of the wire and the screening length 
should be important generic quantities for this issue. In 
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FIG. 1: The general setup contains a molecular wire modeled 
as a cylinder of length L and diameter a between two perfect 
conductors with potentials $i and $2- 



this paper we present a simple model calculation that 
takes the relative magnitudes of these variables explic- 
itly into account. We describe the model in Section II, 



present the calculation in Section III and some results 
and discussion in Section [V. 



II. THE MODEL 

The molecular wire is modeled as a cylinder of length 
L and diameter of order er (the exact way in which a en- 
ters into the model calculation is explained below), per- 
pendicular to and connecting between two planar metal 
electrode surfaces. As depicted in Fig. [jj the cylinder is 
oriented parallel to the z axis, with its axis going through 
the origin in the xy plane. The two electrodes are as- 
sumed to be ideal conductors implying a constant poten- 
tial on the entire surface of each electrode. We set the 
potentials at the left and right wire-electrode interface 
to be <&i = A/2 and $2 = —A/2, respectively. In view 
of (|^) below, this guarantees a vanishing mean potential 
in z-direction and thus a neutral molecule. Finally, we 
restrict the discussion of the potential profile to block- 
ing junctions between electrodes and molecule so that no 
current is flowing. 

As in Ref. || we assume that the wire material does 
have screening capacity, and is characterized by a screen- 
ing length A. It should be noted that the existence, 
nature and size of such screening length in molecular 
systems is an open question with probably non-unique 
answer. Molecules with large gaps between their high- 
est occupied and lowest unoccupied molecular orbitals 
(HOMO and LUMO) will obviously screen poorly, while 
highly conjugated molecules with low HOMO-LUMO gap 
will screen relatively well. 

In the present discussion we assume that over the rel- 
evant length scales (of order L) screening is described by 
a Poisson equation 



According to the cylinder symmetry of the molecule, the 
charge density p(r\\ , z) depends on the radial distance r\\ 
from the wire axis and the position z along the wire. In 
transversal direction, the charge density is assumed to 
be determined by a given molecular electron distribution 
represented by a factor F(r\\). The longitudinal part de- 
pends on the potential along the molecular axis. The 
screening is then described by 
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-F(r||)0>(0,z) 



(2) 



which together with (Q) will allow us to determine the 
potential profile. 

Any assumption about the functional form of F(r\\) 
is in fact an assumption about the confinement of the 
molecular charge distribution in the molecular cylinder 
and in our generic model it is sufficient to take a form 
that reflects the molecular thickness a. Other details of 
F(r\\) are expected to be of secondary importance. 

In the three-dimensional Thomas Fermi model for 
screening in a gas of electrons with charge e and mass 
m e , the screening length A of @ is related to the elec- 
tron density n by 



E F 
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with the Fermi energy 



E F 



(3ir 2 n) 2 / 3 k 2 



2m r 



(3) 



(4) 



At metallic electron densities A is typically of the order 
of 1 A. To have efficient screening in a molecular system 
electrons (or holes) must occupy molecular states that 
are effectively delocalized over the length of the bridge. 
Charge doping by transfer from the metal electrode to 
the molecular bridge may be one source of such elec- 
trons. Their density is expected to be considerably lower 
than metallic, implying a larger characteristic screening 
length. We expect that a calculation based on (|l|) and 
(g) that uses metallic electron density to estimate A will 
provide an upper bound on the effective screening in a 
molecular wire. 



III. THE POTENTIAL DISTRIBUTION 

Using the model described in the previous section, our 
problem is to solve the equation 



V 2 $(r||,z) = ^(r||)$(0,z) 



(5) 



in the range < z < L subject to the boundary condi- 
tions 



V 2 $ = -47T/9 . 



(1) 



$(rii , 0) = A/2 and <J>(ri, , L) — —A/2 . (6) 
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It is convenient to decompose the full potential 

$(r lb z) = $ (z) + (t>(r h z) (7) 
into a first term describing the bare potential 



*o(*) 



(8) 



in the absence of a molecule and a second term which 
reflects the additional potential 4>{r\\ , z) satisfying the 
boundary conditions 4>(r\\ , 0) = <f>(r« , L) = 0. 
The resulting differential equation 



V 2 0(r||,z) = ^F(r||)[$ o (z) + ( /)(O,z)] 

may be solved by the Fourier ansatz 
d 2 fcii 



4k II -r 



(2tt) 2 



^2 $ n (k\\)sia 



(9) 



(10) 



After expressing the bare potential profile in terms of a 
Fourier series one arrives at 



1 



X 2 [k 2 



A 



! + (-!)" 



F^ll) 



(2^ 



.(Mi) 



(ii) 



where 



d 2 r||e- ik " 



\F(r 



(12) 



For the potential profile along the molecular axis, only 
the transversal integral over the Fourier coefficients 



</>„ (fc|| ) is needed which may easily be obtained from (11). 
Due to the symmetry of the bare potential only even 
Fourier coefficients are found to contribute. We thus ar- 
rive at our main result describing the potential profile 
along the molecule 



$(0,2) = $o(z) 
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Jn(l + F„) V L 



(13) 



The coefficients F n accounting for the influence of screen- 
ing are given by 



F n = 



1 



d 2 fc 



X 2 J (2n) 2 k 2 + (2nn/L) 2 



i r, . , /2t™ 

= — / dr||r||F(r||)X (— r„ 



A 2 



(14) 



where Ko denotes a modified Bessel function. In the 
limit of very small screening length, A — > 0, it is possible 
to show by evaluating the sum in ( |l3| ) that the potential 
along the wire vanishes and the entire voltage drop occurs 
at the interface with the electrodes. 



For the following discussion, it is convenient to intro- 
duce a measure of the deviation of the voltage profile 
&(0,z) from the linear behavior (||). Since the integral 
over <&(0,z) — &o{z) vanishes for a neutral molecule, we 
use instead 
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d*($(0,z)-$ o (z)) S 



1/2 



(15) 



This quantity is normalized such that it equals 1 if the 
voltage drop occurs entirely at the ends of the molecule 
while it vanishes for a linear potential profile. Employing 
( |l3| ) , one may express S in terms of the coefficients defined 
by (0) as 



5 = 



6 l/2 



F 2 



"I 1/2 



(16) 



IV. RESULTS AND DISCUSSION 

We now address the dependence of the potential pro- 
file on the width of the molecular wire and start with 
the limiting case of an infinitely thick molecule or, equiv- 
alently, a large number of molecules in parallel present 
between the two electrodes. Then, F(r\\) = 1 and one 
finds from (Hj) 
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Using 



E 



sin 



(nx) it sinh (a(x — 7r)) 



-^n(n 2 +a 2 ) 2a 2 sinh(a7r) 

( |l3| ) yields for the potential profile 

f L - 2z\ 



2a 2 



(17) 



(18) 



sinh 



*(0,*) 



2 sinh(L/2A) 



(19) 



The deviation from the linear voltage drop can be quan- 
tified by inserting ([Tt]) into ([Hi]). Evaluating the sum, 
one finds 



1 



1 + 24 A_- 9 A coth Q^ 2sinh(i/2Aj 



1/2 



(20) 



This result is shown in Fig. |2| as uppermost curve. In the 
limit of very large screening length, A — > oo, S vanishes, 
thereby indicating the expected linear voltage drop. On 
the other hand, for very short screening length, A — > 
0, 5 approaches one and the entire voltage drops at the 
interfaces between wire and electrode. 
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FIG. 2: The deviation S (cf. @) of the voltage profile from 
the linear behavior is shown as a function of the ratio of wire 
length L and screening length A. The four curves correspond 
to a -» oo (cf. @), a/L = 0.05,0.02, and 0.01 from the 
upper to the lower curve. 
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FIG. 3: The potential profile along a molecular wire computed 
from ( 13|) and ( p^ ) is shown for a screening length X/L = 0.05. 
The thickness parameter er/L takes the values 0.0125,0.05, 
and 0.5 from the almost linear behavior to the voltage profile 
containing almost a plateau. 



For the case of finite width, we employ the Gaussian 
charge distribution 



nr||)=exp(-rf/2a 2 ) 



(21) 



This function is not normalized and therefore describes 
a charge distribution with a density in the center inde- 
pendent of the width a. Such a situation arises when the 
diameter of a molecular layer can be controlled. Then, 
the charge density in the center appears in the screen- 
ing length A. In contrast, in the somewhat unrealistic 
case where the charge density on the wire is changed, the 
function F{r\\) would have to be normalized. 

One advantage of the Gaussian distribution ( pl| ) is the 
fact that the coefficients ([t4]) may still be expressed an- 
alytically in terms of an exponential integral 
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with 



1 ( n ° 
X = 2 (2 ? m I 



(22) 



(23) 



With this result the potential profile can be evaluated 
numerically according to (|l^) while the deviation from 
the linear voltage drop is obtained from (|l6|). 

In Fig. [| the deviation S of the voltage profile from 
the linear behavior (^) is shown for different values of 
the wire thickness a. The uppermost curve corresponds 
to the limit of a thick molecular layer a — > oo which was 
discussed above. The three other curves correspond to 
a/L = 0.05,0.02, and 0.01 from top to bottom. As these 
results demonstrate, a reduction of a causes a reduction 
of S indicating that the voltage profile approaches the 
linear voltage drop. This behavior can be understood 
in terms of a reduction of screening due to the reduced 



molecular layer. However, this effect becomes only rel- 
evant for a < I. As discussed above, the limit A — » 
leads to a constant potential along the molecular wire. 
Therefore, all curves shown in Fig. |^ tend to § = 1 in 
this limit even though this is not shown in the figure. 

We now turn to a discussion of the voltage profiles 
themselves. Figure || depicts results obtained from ( fl3| ) 
using F n from ( p2| ) and $o(z) as defined in (||). The 
dimensionless screening length X/L = 0.05 implies for a 
typical metallic screening length A = 2 a.u. a wire length 
of L = 20 a.u.. The thickness parameter a/L for the 
three different curves are 0.0125, 0.05, and 0.5, where the 
voltage profile becomes more and more linear as a de- 
creases. As already mentioned, this may be understood 
in terms of the reduced screening. Fig. [l| shows similar 
results for a wire with a ratio between the typical di- 
ameter and the wire length of a/L = 0.125. Here, the 
dimensionless screening length takes the decreasing val- 
ues X/L = 0.25,0.1, and 0.05 with increasing deviation 
from the linear voltage profile. 

In Fig. H we present a fit of our voltage profile to the ab 
initio calculation of Ref. 0. A least square fit resulted 
in a screening length of X/L = 0.052 and a wire width of 
a/L = 0.032. For a distance L of approximately 34 a.u., 
this yields the reasonable value of 1.1 a.u. for the radial 
extent of the electron density. 

The parameters a/L and X/L employed allow for a 
rather good approximation of the results of the ab ini- 
tio calculation. However, there are two significant dif- 
ferences. The Friedel oscillations found in the quan- 
tum calculation cannot be obtained within our classical 
approach. In a tight-binding model description of the 
molecular wire, Friedel oscillations naturally arise from 
a breaking of electron-hole symmetry jllj . Secondly, the 
systematic shift between the two voltage profiles in Fig. || 
indicates that the wire in the ab initio calculation was 
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FIG. 4: Same as Fig. ^ but for fixed thickness parameter 
a/L = 0.125 and varying screening length A/L — 0.25,0.1, 
and 0.05 from the almost linear behavior to the voltage profile 
containing almost a plateau. 



walled carbon nanotubes of thickness ~ 20 a.u. is con- 
sistent with the results shown in Fig. ||. One should keep 
in mind however that apart from its intrinsic simplicity, 
the model used in this work suffers from two important 
shortcomings. Firstly, the use of a simple screening prop- 
erty as described by ([!]) and (g) cannot be justified for 
all molecules, and certainly not for arbitrary distances. 
Even when such screening applies, the magnitude of the 
screening parameter A is not known and is expected to 
depend on the amount of equilibrium charge transfer be- 
tween the wire and the metal leads. Secondly, a complete 
calculation of the potential profile along a molecular junc- 
tion should take into account the fact that some of this 
drop may take place on the metal leads near the junction. 
Such a behavior was found in the calculation of Ref . la] . 



V. CONCLUSIONS 
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FIG. 5: The potential drop across a chain of 6 gold atoms 
placed between two gold electrodes as obtained from an ab 
inito calculation in Ref. q (dashed line) is fitted by a voltage 
profile (|l^) depicted as full line. The best fit was obtained for 
a/L = 0.032 and A/L = 0.052. 

charged while our wire is always assumed to be neutral. 

It is obvious from these results that the ratio of the 
wire thickness to the screening length does constitute an 
important generic attribute that determines the general 
behavior of the potential bias distribution along a molec- 
ular wire. The relatively good fit obtained between the 
model calculations and the ab initio results for a chain of 
gold atoms using reasonable geometric parameters sup- 
ports this conclusion. We note in passing that the flat 
potential distribution observed pi for metallic single 



The potential distribution along molecular wires in bi- 
ased molecular junctions is determined in principle by 
the detailed electronic structure of the wire and by the 
response of this structure to the molecule-lead contacts 
and to the bias. The present study has identified the wire 
thickness as one of two generic attributes that largely 
determine the way the potential drops along the wire. 
Increasing this parameter leads to a crossover from a 
three-dimensional electrostatic problem to an effectively 
one-dimensional situation. The accompanying increase 
in screening causes a transition from a linear potential 
profile to a situation where the potential drops mostly 
at the interfaces between wire and electrode. The other, 
less accessible molecular property is its ability to screen 
a local charge. In the present model calculation we have 
used a simple screening length parameter to model this 
molecular property, but further studies are needed for a 
better characterization of this property. 
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